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O ABSTRACT 
O 
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^ ■ We report the first results from deep ACS imaging of ten classical globular 

Q ■ clusters in the far outer regions (15 < R p < 100 kpc) of M31. Eight of the clus- 

ters, including two of the most remote M31 globular clusters presently known, 
are described for the first time. Our F606W, F814W colour-magnitude diagrams 
extend ~ 3 magnitudes below the horizontal branch and clearly demonstrate that 
• the majority of these objects are old (> 10 Gyr), metal-poor clusters. Five have 

[Fe/H] ~ —2.1, while an additional four have —1.9 < [Fe/H] < —1.5. The re- 
(Sl ■ maining object is more metal-rich, with [Fe/H] ~ —0.70. Several clusters exhibit 



the second parameter effect. Using aperture photometry, we estimate integrated 
luminosities and structural parameters for all clusters. Many, including all four 
clusters with projected radii greater than 45 kpc, are compact and very lumi- 
nous, with —8.9 < My < —8.3. These four outermost clusters are thus quite 



Or 



2 ■ unlike their Milky Way counterparts, which are typically diffuse, sub-luminous 

(-6.0 < M v < -4.7) and more metal-rich (-1.8 < [Fe/H] < -1.3). 
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Subject headings: galaxies: individual (M31) — galaxies: halos — globular clus- 
ters: general 



1. Introduction 



Globular clusters observed in M31 provide the closest example of a globular cluster 
system belonging to a large external galaxy. Members of this system allow us to trace the star 
formation history, chemical evolution, kinematics, and mass distribution in different regions 
of M31, and are therefore vital to our developing picture of its formation and evolution. 
However, very few clusters have been discovered at large distances from M31, meaning 
that up until now we do not possess a good sample with which to probe its outer halo. 
Rec ently, as part of a ma jor survey of the M31 halo with the Isaac Newton Telescope (see 



e.g. 



Ferguson et al.l 120021 ) and MegaCam on the Canada- France-Hawaii Telescope (see e.g., 



Martin et al.lboOot l. a search for previously unknown globular clusters, at large projected radii 
{R p ) from the center of M31, has been conducted. This search has result ed in the discovery 
of a significant number of new cluste rs, with R p between ~ 15 — 120 kpc (IHuxor et al.l 12004 . 
20051 : lHuxorll2006l ; iMartin et al.ll2006l ) . A sample of 14 of these objects has been the subject of 
deep follow-up imaging using the Advanced Camera for Surveys (ACS) on board the Hubble 
Space Telescope (HS T). Preliminary res ults for four members of the pop ulation of luminous , 



extended clusters of IHuxor et al.l (120051 ) have recently been described (IMackey et al.l 12006 



hereafter Paper I). In this Letter, we present results for the remaining ten clusters, and 
compare, for the first time, the cluster population of the outer M31 halo to that of the outer 
Milky Way halo. 

This sample consists of classical (compact) globular clusters with projected radii dis- 
tributed in the range 15 < R p < 100 kpc. Only two of these objects are previously recorded - 
the remaining eight are new discoveries. Their positions are listed in Table [TJ and presented 
schematically in Fig. [TJ Two of the newly-discovered clusters (GC5, GC10) lie at very large 
distances from M31: R p ~ 78 and 100 kpc, respectively. These, along with the new cluster 

6x3), 



of Martin et al 



which has R p ~ 116 kpc, are by far the most remote members of 
M31's globular cluster system found to d ate. Of the two prev iously reported clusters in our 
sa mple, one (GC4) wa s first described by IHuxor et al.l (120041 ) but discovered independently 
by Galleti et aI. J 2005h. who labelled it B514. This object has additional separate ACS imag- 
ing (jGalleti et al.l 120061 ) . The second was observed in the same field as o ne of our luminous, 
extended clusters, and appears in the catalogue of iBarmby et al. (boooh (cluster 298-021). 
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2. Observations and data reduction 

Our observations were obtained with the ACS Wide Field Channel (WFC) under HST 
program GO 10394 (P.I. Tanvir), at various intervals over the period 2005 May 27 - 2005 
September 17. Targets were imaged three times in the F606W filter and four times in the 
F814W filter, with small dithers between subsequent images. Typical total integration times 
were 1800s in F606W and 3000s in F814W. Each cluster was placed at the centre of either 
chip 1 or chip 2 on the WFC, to avoid the inter-chip gap. Drizzled F606W images of two 
representative clusters are displayed in Fig. [TJ 

Details of our photomet ric reduction procedure are provided in Paper I. Briefly, we 



used the dolphot software (lDolphinll2000l ). specifically the ACS module, to conduct PSF- 



fitting photometry o n all images. Output measurements are on the VEGAMAG scale of 



Sirianni et al.l (120051 ) . We used the quality information provided by dolphot to clean 
the resulting detection lists, selecting only stellar detections, with valid photometry on all 
input images, global sharpness parameter between —0.3 and 0.3 in each filter, and crowding 
parameter less than 0.25 in each filter. Field stars are present in all images; however, their 
densities are typically sparse since the clusters are at such large projected radii from M31. 
To eliminate field star contamination from our cluster colour-magnitude diagrams (CMDs), 
we imposed limiting radii from the cluster centers. Because of the sparse fields, these radii 
were typically quite generous (~ 25"). However, two clusters (GC6, GC7) are set against 
heavier background fields and hence required smaller limiting radii of 10". From the resulting 
CMDs (see below), it is evident that field star contamination is not a significant issue - in all 
cases the prinicpal cluster sequences are very clearly visible. Therefore, a more sophisticated 
statistical subtraction is not necessary at this stage. 



3. Analysis 

Fig. [2] shows CMDs for our ten clusters. The photometry reaches ~ 3 mag below the 
level of the horizontal branch (HB), to a limiting magnitude of mF606W ~ 28. All clusters 
exhibit narrow red-giant branches (RGBs) and clearly delineated HBs. Nine of the clusters 
possess rather steep RGBs, indicating they are metal-poor objects. Many of these clusters 
also feature HBs extending to the blue, including broadened regions with colours in the range 
0.1 < "iF606W — ^F8i4W ^ 0.5 which are suggestive of RR Lyrae stars imaged at random 
phase. Such clusters are old, with ages > 10 Gyr. GC7 has a noticeably different CMD 
from those of the other clusters, with a sharply bending RGB and very stubby red HB, 
characteristics indicative of a more metal-rich object. 
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We obtained photometric metallicity estimates using a procedure we previously de- 
scribed and verified in Paper I. Because our sample covers a very wide spatial area and 
range in R p , we could not assume a uniform distance and reddening. Instead, we determined 
the best-fitting combination of [Fe/H], (m — M) , and EiB — V) for a given clu ster by regis- 



tering the ACS/WFC F606W and F814W Galactic globular cluster fiducials of iBrown et al. 



(120051 ) to the appropriate CMD using the F606W level of the HB and the colour of the 
RGB at the HB level. A fiducial of the correct [Fe/H] closely traces the upper RGB on 
the CMD when registered in this way, while a fiducial of incorrect [Fe/H] deviates on the 
upper RGB. The cluster metallicity is estimated by bounding the RGB with a more metal- 
poor and a more metal-ri ch fiducial. We assume a priori that the main body of M31 has 



(m — M)q = 24.47 ±0.07 (jMcConnachie et al.ll2005l ). For each cluster we registered fiducials 



using, incrementally, a range ±0.5 mag about this value (i.e., a system depth of ~ 360 kpc). 
At each (m — M) we used Brown et al's transformations to solve for the E(B — V) which 
aligned the fiducial and CMD HB levels. Applying this reddening, we determined the offset 
between the color of the fiducial RGB at HB level and that of the cluster CMD. The best 
fitting combination of (m — M) Q and E(B — V) for a given fiducial was that which minimized 
this offset. To determine uncertainties in these values we randomly selected ten thousand 
Gaussian deviates about the two HB levels and the RGB colour at HB level (widths were 
typically ±0.1 mag and ±0.01 mag, respectively) and calculated new E(B — V) and (m— M)o 
for each. We derived standard random errors from the resulting distributions. We note that 
there are possibly additional systematic errors present for our measured E(B — V) values, 
due to differences in age between the M31 clusters and the fiducial clusters, and the fact that 
the colour of the RGB at the HB level is mildly age-sensitive. For the majority of clusters 
(with the exceptions discussed below), these errors would amount at most to ~ 0.02 mag 
more than quoted in Table [IJ corresponding to the interval 10 — 15 Gyr. 

Fig. E] shows the fiducial registration, while numerical results are li sted in Table [H Th e 



derived metallicities reflect those adopted for the reference clusters by IBrown et al.l (120051 ). 
Our typical measurement uncertainties are ±0.15 dex. 

The majority of our clusters are confirmed as metal-poor objects. Five have RGBs well 
traced by the M92 fiducial ([Fe/H] = -2.14), and two have RGBs bounded by M92 and NGC 
6752 ([Fe/H] = -1.54). Two more have RGBs matched by the NGC 6752 fiducial. The re- 
maining cluster, GC7, is more metal-rich, with its RGB matched by 47 Tuc ( [Fe/H] = —0.70). 
Two clusters have previous metallicity measurements. Galleti et al. ( 20061 ) d erived, photo 



metric ally, [Fe/H] ~ —1.8 for GC4. GC6 has a spectroscopic measurement by lPerrett et al. 



( 120021 ). who found [Fe/H] = — 2.07±0.11. Both values are in reasonable agreement with our 



estimates. 
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W e also compared our derived E(B — V) with those from the maps of ISchlegel et al. 
(119981 . SFD) (listed in Table [I]) . Agreement is adequate, although we note a tendency to 



derive larger E(B — V) than SFD by a few hundredths of a mag - we noticed a similar effect 
in Paper I. Adopting the SFD values forces our (m — M) to be greater by ~ 0.025 mag 
per 0.01 mag difference in E(B — V); however the fiducial registration is often noticeably 
inferior. The discrepancy may be due to spatial reddening variations on smaller scales than 
resolved by the SFD maps (6fl), a systematic error introduced by slightly different ages 
for the M31 and te mplate cluster s (as n oted above), or to a systematic error in the colour 
excesses adopted by iBrown et al.l (120051 ) for their reference clusters. 



Most of our derived distance moduli lie close to the canonical value for M31: (m — M) = 
24.47 (~ 780 kpc). However, our measurements for GC7 and GC9 suggest these clusters may 
be closer by ~ 115 kpc and ~ 85 kpc, respectively. This would render GC7, in particular, 
an unusual object, given its [Fe/H]. In the Milky Way, there are two prominent globular 
clusters at unusually large radii for their [Fe/H]: Ter 7 and Pal 12. Both these clusters are 
only about 70% the age of th e oldest Galactic globulars. If GC7 is similarly young, then 
fitting the IBrown et al.l (120051 ) fiducials is inappropriate. We note that the CMD of GC7 
shows an overdensity of points around m F606W ~ 27.5 and mp606W — m F8uw ~ 0.5 which is 
not present in the other cluster CMDs. The presence of the main-sequence turn-off in this 
region, blurred by observational errors, could explain this bulge. If this is the case, then GC7 
may be as young as ~ 4 Gyr, based on CMDs for LMC and SMC clusters. To check our 
measured E(B — V) and (m — M)n fo r GC7, we transformed our ph otometry to Johnson- 
Cousi ns V and I (ISirianni et al.l 120051) and matched the fiducials of ISarajedini &: Lay den 
( 1997 ) and Rosenberg et al. ( 19981 ) for Ter 7 and Pal 12, respectively. These objects have 
comparable [Fe/H] to GC7. The fiducial fits result in E(B-V) ~ 0.07 and (m-M) ~ 24.3, 
much more similar to the quantities observed for our other M31 clusters. We note that both 
Ter 7 and Pal 12 are unambiguously associated with the accreted Sagittarius dSph. GC7 
may represent a similar scenario in M31 - verification of our measurements would therefore 
be valuable. 

Inspection of Fig. [2] and our measured [Fe/H] reveals several second parameter clusters. 
We defer numerical calculation of HB morphologies until a future paper when tests for 
completeness and photometric blends can be incorporated. Even so, the second parameter 
effect is clear for GC5, which has a red HB but equivalent [Fe/H] to GC2, and EC1 from 
Paper I. Similarly, GC8 and GC9 have [Fe/H] ~ —1.54; however they have much redder 
HBs than comparable Galactic globulars (e.g., NGC 1904, 6752, 7492). The HBs of these 
two objects resemble those for many of the M31 clusters with similar [Fe/H] but smaller R p 
observed by iRich et al.l (120051 ). and second parameter Galactic globulars such as Pal 14. 
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Finally, we calculated integrated cluster luminosities by means of aperture photometry 
centred on the cluster centres. This technique ensures all light is counted, even in the 
unresolved cores of the most compact objects. For a given cluster, we first estimated the sky 
level using regions away from the cluster, and produced a sky subtracted image. At large 
projected radii from M31, we do not have to worry about a significantly spatially variable 
background. We then masked any bright background galaxies or foreground stars in the 
vicinity of the cluster. Next, we measured integrated luminosities using a variety of apertures. 
In combination with the cluster CMDs, plotting luminosity as a function of aperture radius 
allowed us to check for field star contamination. For the most remote clusters, where field 
contamination is extremely sparse, the integrated luminosity quickly reached an asymptotic 
limit. For such objects a maximum radius of 20" was more than sufficient. For objects set 
against more significant fields, contamination manifested in the form of a slowly increasing 
integrated luminosity with increasing aperture radius (as well as being visible in the CMD). 
For these objects, the cluster light always dominates for radii smaller than 10". Beyond 
this, we used our plots to estimate the radius at which field star contamination became non- 
negligible (r max ) and integrated only to this limit. Three clusters have r max = 10", while the 
remainder have 16" < r max < 20". We used our summed magnit udes in F606W and F814W 



to transform the F606W value to Johnson V ( iSirianni et al.ll2005l ). and then used the relevant 



E(B — V) and (m — M)o to calculate My. We next re-sampled with smaller apertures to 
estimate the half-light radius r^, which we converted to parsecs using our (m — M)o values. 
All results are listed in Table [TJ Clusters with r max > 16" have reliable My - as noted above, 
the integrations are very close to asymptotic by such radii. For the three with r max = 10", 
we used the profiles of the clusters with negligible field contamination to observe that My 
may be under-estimated by ~ 0.1 mag. This level of uncertainty is acceptable for present 
purposes - it is less even than the uncertainties introduced by application of the derived 
distance moduli. 

With our new sample of remote M31 members, we are, for the first time, in a position 
to compare the outer globular cluster system of this galaxy with that of our own. Doing so 
reveals striking differences between the two. Seven Galactic globular clusters have R gc > 40 
kpc (Pyxis; Pal 3, 4, and 14; AM-1; Eridanus; and NGC 2419). Of these, all except NGC 
2419 have —1.8 < [Fe/H] < —1.3, and are sub-luminous (—6.0 < My < —4.7) and very 
extended (11 < < 25 pc). NGC 2419 is one of the most luminous clusters in our Galaxy 
(My ~ —9.6), and lies in an unusual position on the luminosity vs. size plane. This has led to 



some suggestions it may be the remains of an accreted galaxy (e.g.. Ivan den Bergh fc Mackey 



20041 ) . In our present sample, we have four M31 globulars at R p > 40 kpc. In contrast 
with the majority of outer Milky Way clusters, these objects are very metal-poor (—2.2 < 
[Fe/H] < -1.8), compact (4 < r h < 7 pc) and very bright (-8.9 < My < -8.3). Unlike 
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NGC 2419, none of the clusters lie in an unusual position on the lumi nosity vs. size plane 



suggesting that they are 'normal' globular clusters. The new cluster of iMartin et al.l (120061 ). 
at Rp ~ 116 kpc is also compact and luminous {My ~ —8.5), while EC4 from Paper I, at 
R p ~ 60 kpc, is a member of the population of luminous, extended M31 globular clusters 
which has no Galactic analogue. 

It is perhaps not surprising that we have not observed any very extended, sub-luminous 
globular clusters in the remote M31 halo - such objects may well lie below our present survey 
detection limit. Even so, the above comparison shows that M31 clearly possesses an extended 
system of metal-poor, compact, and very luminous globular clusters which is not seen in the 
Milky Way. The one luminous, metal-poor outer Milky Way cluster - NGC 2419 - is quite 
unlike the M31 clusters we have observed (as described above). Furthermore, the outermost 
M31 globular clusters studied here are consider ably more metal-poor than the pressure - 



supported field halo population at the same radii (IChapman et al.ll2006l ; iKalirai et al.ll2006l ). 
again in contrast to the Milky Way. These disparities may well offer important clues to 
differences in the early formation and evolution of the two galaxies or in their subsequent 
accretion histories, and as such, it is vital that additional M31 members are sought. We 
remark that the present sample has been obtained from a survey of ~ 100 square degrees, 
roughly 25% of the area contained within a projected radius of 150 kpc of M31. A significant 
number of new clusters may therefore await detection, suggesting that, unlike in the Milky 
Way, there may be a relatively large population of luminous, compact globular clusters in the 
outer M31 halo. Such objects, in addition to their value as globular clusters, are extremely 
useful as dynamical probes of the mass distribution in M31 - work which has previously 
relied on clusters with R p < 25 kpc and a handful of dwarf galaxies. We anticipate that our 
new sample will soon allow improved mass constraints to be determined. 
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Table 1. Observed properties of ten globular clusters in the halo of M31 
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Barmbv et all |2000h . 
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Fig. 1. — The upper panel shows the location of our present cluster sample (open diamonds) 
in relation to M31. The four luminous, exte nded clusters from Paper I are marked with 
stars. Also displayed is the remote cluster of Martin et al.l (120061 ) (open triangle) and two 
M31 companion galaxies - M32 and NGC 205 (filled circles). The small ellipse delineates the 
visual extent of the M31 disk. The lower panel shows drizzled ACS/WFC F606W images of 
the two outermost clusters in our sample (R p ~ 78 and 100 kpc, respectively). Visually, they 
are representative of the classical globular clusters considered in this paper. Each thumbnail 
has dimensions of 25" x 25". 
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Fig. 2. — CMDs for the ten classical M31 globular clusters. Photometry has been selected 
from the full ACS/WFC fields by imposing radial limits as described in the text. Typical 
photometric errors from dolphot are indicated. 
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Fig. 3. — Results of fitting the Galactic globular cluster fiducials from lBrown et al.l (120051 ) to 
the observed CMDs. In all panels, the more metal-poor fiducial is marked with a solid line, 
and the more metal-rich fiducial with a dashed line. The numbers in the upper left of each 
panel indicate the best fitting (m — M) and E(B — V) for the two fiducials - metal-poor 
above, and metal-rich below. The two marked fiducials for GC7 are 47 Tuc and NGC 5927, 
which have [Fe/H] = —0.70 and —0.37, respectively. For GC8 and GC9 the two fiducials 
are NGC 6752 ([Fe/H] = —1.54) and 47 Tuc, while for the remaining seven clusters the two 
fiducials are M92 ([Fe /H] = —2.14) and N GC 6752. Photometry has been transformed to 
STMAG to match the iBrown et al.l (120051 ) measurements. 



